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During the first order phase transitions the objects of a new phase consume 
a mother phase which leads to some density gaps near existing embryos of a 
new phase. This fact is rather obvious [Q, but complete theoretical description 
was given only in ph. Now this fact in widely recognized |B[ but still there 
are some attempts pj to reconsider this approach which is now known as the 
model of depletion zones || . Here we shall analyze Q and show the useless of 
reconsideration proposed in Q . 

We don't consider here the picture of the gobal evolution of a system - in j| 
it is so primitive that can not give any reliable quantitative results but only some 
very rough estimates without any justifications. We needn't to do it because this 
has been already done in [Q] with all necessary justifications. So, wc consider 
only the profile around the solitary embryo, which was the matter of interest 
in Q . Why it is possible to consider only the solitary embryo? An answer has 
been given in in Q| this question remains unsolved. Due to results of |2j we 
shall consider namely the solitary embryo. 

We assume that the substance exchange regime is a diffusion one. Then 
there is a gap of a vapor density n near an embryo. Then the density n differs 
from a density n(oo). This difference initiates a variation in a nucleation rate 
/ and causes the effect of depletion zones. So, it is very important ro get for 
n(r) ~ n(oo) (r is a distance from a center of an embryo) a correct expression. 
This expression is a crucial point in consideration of an effect of depletion zones 
in nucleation kinetics. Certainly, it has to be a function of time t. 

In Q the following expression is proposed 

n(r, t) = n s (r) + (n(oo) + n(0))-*(^=) (1) 

' zv L>t 

, \ , s n(co) — n(0) „ 

n s (r)=n(oo) — —R (2) 

r 

Here n(0) is equilibrium concentration on the surface of nucleus, D is a diffusion 
coefficient, R is the nucleus radius, $(a:) is the Laplace probability integral. 
Here the sign + was used instead of + in the second term of (El). Correct 
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solution in the neigbourhood of an embryo is Q 

n - ra(oo) R r - R 2 f°° 2 

— L' = -er/c — -= , erfc(z) = -= / exp ^ (3) 

But even with this solution one can not construct the correct form of the 
nucleation rate profile. We shall outline this point. 

There are following necessary suppositions to get this exact solution for 
concentration profile: 

• The radius of an embryo is constant in time 

• The boundary condition at r = R is constant in time. 

Both assumptions are rather approximate. As a result the above solution is 
valid only near the surface of an embryo. It is easy to see from the following 
arguments. Let St be the characteristic time when the relative variation of 
intensity of vapor consumption v is already essential 

\v(t + St) - v{t)\ ~ v(t) 

Then it follows that only the distances 

r -R~ ADSt = 5R 

or smaller ones can be considered on the base of solution (||). At all distances 
r — R > SR we have to take into account that earlier the intensity of vapor 
consumption was smaller. Really 

v ~ du/dt 

where v - a number of molecules inside the embryo. As it is known 

dv/dt = const j/ 1/3 - t 1/2 

and we see that the intensity isn't constant in time. 

It seems that the small distances are the main ones. Certainly the gap is 
greater the smaller the distance is. But the small distances aren't essential 
because the functional form of the nucleation rate ^ can be approximated as 

I(r) = I(r = oo) exp(r(n(r) — n(oo))/n(oo)) (4) 

where T is a big parameter. It is approximately equal to v. Note, that the 
sence of nucleation rate as a probability to appear for an embryo was used here 
(details see in So, the distances where 

— r(n(r) — n(oo))/n(oo) ~ 1 

are the most interesting. When 

\T(n(r) - n(oo))/n(oo)| < 1 



2 



we have I — I{oo) and the is no gap of concentration. When 

— r(n(r) — n(oo))/n(oo) 3> 1 

we have 1 = and the is no nucleation. So, this region isn't interesting. 

The distances T(n(r) — n(oo))/ro(oo) ~ 1 because of T ~ v — > oo corresponds 
to very small n(r) — n(oo) and, thus, to big values of distances ~ Sr from the 
embryo. Ordinary Sr ^> SR and there is no possibility to use (|J). 

Instead of ||) in Q another approach was used. This approach is based on 
the strong inequality 

Sr > R (5) 

Why this inequality takes place? When we consider classical nucleation, i.e. 
transition of a supersaturated vapor into a liquid phase everything is clear. We 
have a strong inequality 

Vv/Vl > 1 

where v v is a partial volume for one molecule in a vapor, vi is a partial molecule 
in a embryo. This inequality makes obvious that the final volume after the whole 
process of condensation is small. As an example of condensation of water vapor 
in normal external conditions. Suppose that supersaturation is somewhere ~ 5. 
Then because v v /vi ~ 1000 we see that the final volume of a new phase is very 
small ~ 1/200 and the characteristic distance Sr ~ (200) 1 / 3 i? » R. 

We see that the value Sr has to take the same order of magnitude as the 
mean distance between objects of a new phase r mean . Really, when Sr attains 
imaginary value r mean it means that the process of nucleation stops. So, namely 
the values at this moment are the final ones and the main ones. Then 

Tmean ^° Sv 

The same will be valid for all other first order transitions. 

Consider the opposite transition from the liquid to the vapor phase. We 
shall see the opposite inequality. But we know that compressibility of liquid 
is very law. So, the relative super streching is rather low. It means that the 
vapor phase in a final state will occupy some rather small volume. And again 
we see that the mean distance between the objects of a new phase is many times 
greater than the size of objects. Again we get the same result. 

We see that the following fact takes place: The relative final volume of a 
new phase will be small. 

To see that (||) takes place we needn't to consider the real final volume of a 
new phase after phase transition but can take the values at the end of nucleation. 
Because the relative quantity of surplus substance condensed in a new phase has 
an order T _1 (at least under the collective regime of substance consumption 
we see that the relative volume occupied by a new phase is limited from above 
by a value of an order r . So 

r - i?r 1/3 

' mean Ju 

and the required property is established. 
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So, we see the required property for moderate effects of depleted zones. 
When the effect of depletion zones is strong we can use at r ~ (2~3)R stationary 
solution and see that Sr > (5 -f- 7)R. So, here the required property is also 
observed which completes justification. 

Then we can use profiles obtained on the formalism of Green function. We 
have § 

n(oo) - n(r) = / ^ --^ exp(-— ^ — — )dx (6) 



o 8(Dir(t - x)f/ 2 lK AD{t-x) 

A^(4^) 3 / 2 (^) 1 / 2 (Cn 00 ^) 3 / 2 (7) 
Ztt 

and C is a supersaturation. 

This result was obtained for transition from vapor phase into a liquid phase 
but there is no difference because the kinetic mechanismremains the same. 

One can see that expressions (|l|), (Q) and (^), (R) are absolutely different in 
analytical structure. The same is valid when we compare (Q), (0) and (^), (0) 
numerically. 

On the base of profile we can easily get all results of the nucleation process. 
The rate of nucleation can be found according to ([|). Then we have 

I(r) = I(r = oo) exp(rn(0) — n(oo) — er/c(— ^=)/n(oo)) (8) 

for analytical solution with fixed boundary of embryo and 

J(r)=J(r = oo)exp(-ry o ^ - ^ expC-^— y )dx/n(oo)) (9) 

for solution on the base of Green functions. Note that the analogous transfor- 
mation in p| was done with partial shift of coordinate r — Rtor which is wrong 
and beside this isn't necessary. 

We see that these approaches give absolutely different results. 

The next step is an obvious remark that the probability dp that an embryo 
appears during elementary interval dt at the distance from r up to r + dr from 
the center of the already appeared embryo is 

dp = dt^r 2 I(r,t) 

Then the probability dP that the embryo appears in the layer from R up to r is 

dP = dt [ 4Trr' 2 I(r',t)dr' (10) 
Jr 

This expression differs from M] where the following expression was presented 



dP = dt / ATir 2 I(r\t)dr' 

J r 
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The last expression is wrong and we shall follow @). 

The next step used in Q] is to come to an integral value P which is the 
probability that in a sphere of radius r around already existing embryo there 
will be no appearance of a new embryo. Certainly we can start from ( |l0| ) and 
integrate it. But it is necessarfy to take into account that the zone of depletion 
will grow according to the absolutely precise law. This law was established in 
[pi . But here we suppose^ that one can act in another manner and write 

P = exp(- / / 4irr' 2 I(r',t')dr'dt') (11) 



JO JR. 

The last expression analogous to those which forms the base for further 
constructions inQ Q is rather doubtful. The problem appears because now we 
are coming to characteristic of a whole process of nucleation (earlier there was 
a profile around a solitary object of a new phase). 

For a "solitary subsystem" the derivation of ( |Tl]) is evident: Let 

dp = dtl 

be a probability that during an elementary interval dt there will be appearance 
of a new embryo. The value of dt is rather small. Then dp is proportional to 
dt and to the rate of nucleation I (here the sense of the nucletaion rate as a 
probability is used). Then the probability for the absence of appearance is 

dp' = 1 - dp = 1 - dtl 

If we have two elementary intervals dt\ and dti which don't overlap then the 
probability of the absence of embryo appearance is 

dp\dp' 2 = (1 - 7i<Ki)(l - hdt 2 ) 

The same is valid for arbitrary number of intervals. 

To fulfill multiplication in the r.h.s. one can present dp' as 

dp' = (1 - Idt) = cxp(-Idt) 

Then the total probabiblity of the absence of appearance is 



and we come to (jlj 

The real problem is how to take into account the overlapping of density 
profiles initiated by different embryos. This problem was solved in The 
exponential form (Thy) for P is some rather arbitrary interpretation of the value 



1 To be close to | 
2 In Q an interga 



over r is from up to r. 



5 



"the free volume" used in 0. Now we shall discuss this interpretation in frames 
of the approximation of solitary droplet. 

The overlapping of exhausted regions (ER) requires to use instead of (??) - 
(^another approximation where instead of one profile there is a superposition 
of profiles. Certainly it is impossible to calculate this superposition precisely. 
Then one will come to some models analogous to formulated in . 

One can directly use the results from j2| here. When the free volume Q is 
calculated as function of time one can determine the total nucleation rate as a 
ratio of the free volume to the whole volume of a system and use then (Jl|). 

The real problem is what consequences can we make from a knowledge of P. 
On one hand it seems that the knowledge of P solves al problems in knetics of 
the global nucleation stage. Really, having presented P in a form 

P = exp(-L(r,i)) 

where L(r,t) is some expression one can approximately estimate the average 
radius of a sphere where there will be no appearance by relation 

L(r,t-+ oo)| r=ro = 1 

This construction analogous to needs two remarks. 

This estimate makes no difference between embryos appeared earlier or later 
during the nucleation stage. It means that this approximation supposes that 
the spectrum of droplets sizes is monodisperce one. It isn't correct assumption 
and can be treated only as a rough estimate. 

It is known that sometimes the overlapping of profiles is important and 
makes the main contribution in nucleation kinetics. This phenomena can occur 
when the substance exchange regime is goimg to be a free molecular one and 
when the long tails of profiles are important. The first situation can not be 
realized here. The last situation was a metter of separate consideration in . 

In (|] an integral definition of a boundary of ER allowed to take into account 
the situation of long tails. Here this situation can not be considered, because 
the level type definition is used. It means that the value of rg is determines as 
a value when L attains some level (=1). The definition of level type can not 
take into account long tails and it is more convenient to use approach from Q . 
Here this approach is considered because it is analogous to ||]. 

When the value of ro is determined it is easy to estimate the total average 
number of droplets N as 

(the coefficient can be omitted, this depends on interpretation of ro). 

To fulfill concrete calculations one can use the following approach. Note that 
P can be in any case presented as 

P = exp(- / dt' dr'Anr' 2 exp(/(r',i')) 

JO Jr l 
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The lower limit r/ of integration has to be put n = R for models with explicit 
boundary of an embryo and for models with Green function r; = 0. But because 
5r S> R the effect of the lower limit will be small and we can put n = in 11 
situations. Certainly, for P due to the exponential form the effect is essential, 
but when we are calculating the mean distances the essential dependence on r; 
disappears due to Sr ^> R. In the last equation / is some function with explicit 
form given by solution of diffusion equation. 

One can see that in f(r,t) two variables has to appear in combination (3 = 
t/r 2 or j3' = t/(r — R) 2 . The dependence on t and r via j3 is the main one (in 
special regimes of mother phase consumption there may be dependence on t, r 
in another combination, but this dependence will be much more smooth than 
via fi). Because SR^> R the difference in substitution j3' instead of @ will be 
small (not in P, but on mean distances and times). 

Then in integration J Q dt' we have to substitute dt' by d(3' which gives 

P = cxp(— / dr' 4-Kr' 4 const 
Jo 

where const comes from 

/•oo 

const = / d/?'/3' 2 exp(/(r',i')) 
Jo 



because 



poo POO 

/ dt' = r' 2 / dp' 
Jo Jo 



It leads to 

P = cxp(-r 5 7 ) (13) 

where 7 = 4nconst/5 is some constant. 
One can show that 

/ dxexp(— x 5 ) w 1 
Jo 

This approximate equality is very important. Namely this equality allows to 
determine the mean distance between droplets as TV" 1 / 3 , or as tq or according 
to the sense of P. Three ways are available. They have to give similar results. 
This fact is ensured by this approximate coincidence. Then the third way gives 



/■OO (~1 p 00 



where Co is the normalizing factor of distribution P. One can calculate Co 
according to 



/>oo 

C / drexp(-7r 5 ) = 1 
Jo 



and 

r 7 1/5 



J dxcxp(—x 5 
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Theif] 

1 J °° c?aexp(-a 5 / 2 ) 
27 1 / 5 J °° dx exp(— x 5 ) 

Because 

cfaexp(-a 5 / 2 ) w 1 



/>00 

/ dxexp(— a; 5 ) 1 
Jo 



we see that r is two times smaller than the mean distance between embryos. It 
is clear because the mean distance between embryos is two mean distances until 
the boundary of exhaustion zone tq. So, we observe a coincidence. Certainly, 
one can calculate integrals more precisely. 

The average time of waiting fro the appearance of new embryo near the 
already existing one is directly and elementary connected with a mean distance 
between embryos and we needn't to calculate it separately^. 

As for calculation a quadratic mean deviation from the average distance it 
is given0 by 

< (r- <r>) 2 >=<r > 2 & - 1) 

where 

/■OO 

Ii = I dxx 1 exp(— a; 5 ) 
Jo 

But we have to stress that this result has practically no meaning because as it 
is stated in Jt| the effect of interaction leads to decrease of fluctuation of the 
total number of droplets. The right numerical value of this effect is given in ||. 
Here we have absolutely no interaction between droplets and this result has to 
be seriously changed due to interaction. An example of such account is given in 



To close this question we shall construct now more precise solutions of dif- 
fusion equation. They have to be used in (|]) and then in ( p"l| ) 

We shall start with the following model. There exists a more general solution 

0. 

For a domain R < r < 00 where 

n = f(r) at t = 

and 

n = g(t) at r = R 
the solution of diffusion equation 

dt dr 2 r dr 



3 This value is calculated in M in a wrong way. 

4 The integral in expression for this value is calculated in M in a wrong way 
5 Result in W is wrong. 
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(this the diffusion equation with spherical simmetry) is || 



2r\ n KDt Jr 



where 

r- R 



(14) 



+ ^ / git - V > ) exp(-r 2 )dr 



This solution allows to formulate the following approximation. Let R be 
constant 

R = R{t = 0) 

but the boundary conditions will be recalculated according to the stationary 
solution: 

n s (R(t = 0)) = n(oo) - ^I^W ) (15) 

Really the rate of embryo growth is so small that the stationary distribution 
can be regarded as practically precise one. So, this approximation is practically 
precise and takes into acound the variation of mother phase consumption by an 
embryo. 

Note that diffusion equation 

9 ™=D\ 



dt dr 2 r dr 

in the system with spherical symmetry can be transformed by substitution 

u(r, t) = rw(r, t) 

to an ordinary one dimensional diffusion equation 

du jj® 2u 
dt dr 2 

We shall consider this equation in future. 

Then we can make a shift x = r — R(0). The problem 

w = f{x) 

at t = 

w = g(t) 

at x = has solution 



(16) 



2V^7o V 4Z?(i-r) y (t-r)»/a 
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It is also possible to go to variable 

p = r- R(t) 

where R(t) = At 1 / 2 is a boundary of an embryo. Then in the diffusion equation 
beside 

- _ 1 du d 2 u 

there appear some additional terms. Having considered Q as a main operator 
we come to an iteration procedure where all other terms are assumed to be 
known in a previous approximation (in zero approximation they are omitted). 
So, at every step of iterations we have to solve a problem 



at t = 

at p = R(0) for equation 



to = f{x) 
to = g(i) 



1 du d 2 u 



may be with renormalized D. Here $' is a known (at the previous step) function. 
After the mentioned shift x = p — R(0) we came to a problem 



at t = 

at p = R(Q) for equation 



to = f(x) 
w = g{t) 

1 du d 2 u 



with known The solution of the last problem is [|| 



u= / G(x,£,t)f(Odt 



2^Jo AD{t-rY{t-rf/ 2 



ex P[ — 77T77 ZT 1 7T^~^T : )A T W 



where 



Jo Jo 



The last relation solves the problem to construct iterations on every step. Thus, 
we come to solution of diffusion equation which has to be used in our previous 
constructions. 
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